function K = faster_kernel_function(d1,d2)

ell = 1;                                 
sf2 = 1;                                          
K = distance_gpr(d1'/ell,d2'/ell);
K = sf2*exp(-K/2);

end